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Abstract 

In this paper the attractors synthesis algorithm for a class of dissipative dynamical systems with 
hyperbolic equilibria, presented in [l] , is applied to generate any attractor of the Rikitake system. By 
switching periodically, or even randomly, the control parameter inside a given set of values, during any 
finite time interval while the attractor is numerically approximated, any attractor can be generated. 
Beside the extension of the synthesis algorithm to systems with non-hyperbolic equilibria, we have 
found for the Rikitake system, a new intriguing transient which, occurring for a long time interval, is 
difficult to be numerically found due to the known system instability along the xs-axis. 
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1 Introduction 

The paleomagnetic records of the Earth's magnetic field show that the field has changed its polarity many 
times along geological history (hundreds of times during the last 160 million years). But intervals among 
such geomagnetic polarity reversals are highly irregular. Thus while their average is about 7.10^ years, 
there are intervals as long as 3.10^ years without polarity change, but with large deviations of the poles 
from actual positions. 

There are many factors which can affect the Earth's magnetic field. For example the Reynolds number 
of the Earth's liquid core is believed to be of the order of 10^, i.e. sufficiently large that the flow of electrical 
currents in the liquid core be turbulent. Hence the pattern of the Earth's magnetic field is very complex 
(see for ex.[51|3]). Because of such great complexity, to study the reversals of the Earth's magnetic field, 
relatively simple mechanical dynamos had been proposed as analog models [H [SJ |B] . One of them is the 
Rikitake dynamo of two frictionless coupled disks, which is a paradigm of the geomagnetic field behavior, 
proposed by the Japanese geophysicist Rikitake j4j. This model can be considered only as a special case 
for the real geodynamo, which is obviously a high order physical system with very large degrees of freedom 
[71 [5] . The Rikitake dynamo is composed of two conducting rotating disks which are connected to two 

coils so that the current in each coil feeds the magnetic field of the other (see FigfTl). Each circuit has 



the same self-inductance L and electrical resistance R, and for each disk dynamo a constant mechanical 
torque G is applied from outside on the axis, so that it can rotate with the angular velocity w. 

The currents Ii and I2 in the circuits and the related voltages Vi and V2 are connected by the well 
known relations for the R, L circuits 

Rh+L^ ^ Mhioi, 
Rl2 + L^ ^ MI2UJ2, 

where MIiUJi and MI2UJ2 are the voltages Vi and V2 , and M is the mutual inductance. If C is considered 
the moment of inertia of each disk, these relations can be rescaled and lead finally (see for ex [5]) to the 
model equations 

Xi — X2 * X3 — axi, 

2^2= {x3 - p)xi - ax2, (1) 
2:3= I-X1X2, 

with a = R\J and p = (cji — UJ2) \J The control parameter p is considered to be positive. 

The system is invariant under the change {xi,X2, x^) 1 — > — X2, X3). Therefore if {xi{i), X2{t), X3{t)) 

is a solution to ([T]), then {—xi{t),—X2{t),X3{t)) is also a solution. 

The system is dissipative the divergence being negative divf{x) — — 2p, where the vector function 
/ : ^ jg ^-jjj, j-igj^i- Yiand side of 

The equilibrium points are 

Xi,2 (^±k,±^,ak''^ , with fc' = ^ (p + Vp^ + 4a2) , 

and they are not hyperbolic. 

The Rikitake model is still intensely investigated, not only for its physical interest, but especially for 
the richness of its dynamical behavior (see for example |10| , lll j ). In our numerical research of the dynamics 
of this system, we have found an interesting 'transient attractor' (TA) persisting for a long time interval 
before the trajectory reaches one attractor. The trajectory generating this transient passes through the 
'real' attractor several times. TA size is disproportionately high compared to the size of real attractor. 
Moreover, it is interesting to see that TA is very unstable being related to the xa-axis. 

For that reason our paper is focused on two directions: primarily, to present TA with its proper 
behavior and secondarily to synthesize any Rikitake attractor by using the attractor synthesis algorithm 
(introduced in Jj), that switches the control parameter of the system for finite time intervals, while the 
model is numerically integrated. Therefore, this algorithm firstly applied to systems with hyperbolic 
equilibria, is extended here to another class of systems. One of the main algorithm benefits is the fact 
that it allows the generation of any attractors, because of the convexity property presented in Section 2, 
even if for some objective reasons, some parameter values are not accessible. 

The organization of the paper is as follows: Section 2 describes the synthesis algorithm while in Section 
3 the TA is presented beside the application of the synthesis algorithm. The Conclusion summarizes the 
results of this work. 
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2 Synthesis algorithm 



Let us consider the following Initial Value Problem (IVP) 

x{t) = f{x{t))+p{t)Ax{t), x{0)=xo, ie/=[0,oo), (2) 

where / : M" — ^ M" is a nonlinear vector function, xq E M" and p : I M" is a piece-wise continuous 
periodic function with period T and mean value p* , i.e. 



t+T 

p{u)du = p* , t £ I, (3) 



and A is a real n x n squared matrix. 

Next, we make the following assumption 
HI The IVP ^ admits unique solutions. 

In [12] it is proved that the solutions of the IVP ([2| for the case of continuous systems with respect to 
the state variable, and of the corresponding averaged model, expressed as follows 

xit)^ f{xit))+p*Axit), x{0) = xo, tel, (4) 

may have arbitrarily close solutions. This means that switching p in the IVP ^ following some periodic 
scheme within a selected set of values, while the IVP is numerically integrated, the solutions remain close 
enough to the solutions of the IVP (4). The analysis is carried out based on the averaging theory [T5] . 

The class of systems modeled by (2) includes known dynamical systems such as: Lorenz, Rossler, Chen, 
Lotka-Volterra, Lii, minimal networks, neuronal networks, a class of lasers, etc. Moreover, by using several 
computational tests, we verified that the switching algorithm applies not only to continuous dynamical 
systems, but also to fractional-order system^ 

While in the mentioned examples we have studied systems with hyperbolic equilibrium points, in this 
paper we prove numerically and computationally, that the algorithm can be applied to the Rikitake's 
system which has non-hyperbolic equilibria. 

Despite the fact that could be some differences between computation and theory, the numerical in- 
tegration of ([2| can generally give excellent approximations to the trajectories within the invariant sets. 
The trajectories that start near an attractor will stay near and they will be shadowed by orbits within 
the attractor because attractors arise as the limiting behavior of trajectories. Therefore, the shadowing 
property [16j enables us to recover long time approximation properties of numerical trajectories necessary 
in our numerical computations. 

Roughly speaking, a global attractor can be viewed as a region of a dynamical system's state space 
where some of the system trajectories can enter and not leave, and which contains no smaller such region 
[T7] . A global attractor contains all the dynamics evolving from all possible initial conditions. In other 
words, it contains all the solutions, including the stationary and periodic ones, as well as chaotic attractors, 
relevant to the asymptotic behaviors of the system. 

The term local attractor is sometimes used to denominate non-global attractors (see e.g. [H])- The 
global attractors may contain several local attractors. Therefore, a global attractor may be considered to 
be composed of a set of all local attractors. Each of them only attract trajectories from a subset of initial 
conditions, specified by its basin of attraction. 



^For example, the switching algorithm was applied successfully to a fractional variant of Lii system .14^ . whose chaotic 
behavior is analyzed in I15| . 
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Remark 1 For the sake of simplicity, hut without loss of generality, in this paper when a global attractor 
is composed by several local attractors, only one of the single local attractors will he considered. 

The attractors are numerically approximated using some scheme for ODEs with fixed step size h, after 
the transients are neglected (see e.g. ( i¥) ) . 

We have shown via numerical approach and computer simulations, that while the model ^ is numer- 
ically integrated, if one switch the control parameter at finite time intervals, the obtained approximated 
attractor {synthesized attractor) is approximately identical to the one corresponding to p — p* {averaged 
attractor) for whatever considered set of values for p. The algorithm, we call hereafter synthesis algorithm 
(SA), consists in using a time-periodically parameter switching, according to some designed rule. It will 
be demonstrated empirically, that any desired attractor can be duly obtained by the proposed switching 
scheme. Moreover, we found out numerically and computationally that SA can be applied not only via 
some periodic parameter switching rule (as analytically proved in [12 ), but using any kind of random 
switching way. 

Notation 2 Let A he the set of all attractors depending on the parameter p, including attractive stable 
fixed points, stable limit cycles and chaotic attractors; let V CM. be the set of the admissible values of 
p and Vn = {pitP2t ■ ■ ,Pn} d V a finite ordered subset of V which determines the set of attractors 
An = {^pi7^p2i • ■ ■ i^pn} 

Remark 3 It is natural to introduce a bijection between the set of all admissible values of p and the set 
of all attractor^ Therefore, giving any p, a unique attractor is specified, and vice versa. Also, via this 
bijection, the order over V induces an order over A. 



For the sake of simplicity, unless necessary, we denote the attractors Ap. simply by Ai. 
We assume that we can access all the values of Vn — {pi,P2, ■ ■ ■ ,Pn} for which the system behaves 
stably or chaotically. 

With a chosen finite subset Vn, the SA relies on the following deterministic time switching rule applied 
indefinitely on /, while a numerical method with fixed step size h integrates the IVP ^ 

[pi|/i,P2|/2, • ■ ■ >-PJv|/jv]: Pi ^"Pn, i = 1,2,. ..,iV, (5) 

where li, i G {1, . . . , N} arc finite consecutive (adjoint) time subintervals of length Ati, for i = 1,2, . . . , N. 
([5]) means that in each interval li, p = pi, i — 1,2, . . . , N. In other words, p is a piece-wise continuous 
(constant) function p : li — > Vn, p{t) = Pi, for t ^ li, i = 1,2, . . . , N. 

The simplest way to implement numerically ([5| can be described by the scheme 

[Pi\ mih, P2\m2h, ■ ■ ■ , PNlniNh ], (6) 

where Ati is chosen to have the length Ati — mih with rui positive integers (see Section 3). SA acts as 
follows: in the first time subinterval Ii of length mih, p{t) = pi, then for t € I2, p{t) = P2 and so on until 
the iV-th time subinterval of length mj^h where p{t) = pn. Next, the algorithm repeats. Relation (joj) is 
periodic with T = {mi +m2 + ... + mN)h. In order to simplify the notation, for a fixed step size h, scheme 
^ will be denoted hereafter 

^Thc rigourous proof of this afirmation, remains a future objective. 
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[mipi, m2P2, ■ ■ ■ ,mM Pn]- (7) 
For example, for A'^ = 3, by the scheme [Ipi, 3p2j 2^3] one should understand the infinite sequence of 
p : pi, p2, P2, P2, P3, P3, Pi, ■ ■ ■ which means that while ^ is integrated, p switches in each / subinterval 
between the values of Vs — {pi,p2,P3}. 

Remark 4 In practical examples, switching techniques can be applied not only to parameters but, for 
example, to the state variables 120]. 

In order to compare two attractors, we introduce the following criteria 

Definition 5 Two attractors will be considered approximately identical (AI) if after neglected transients, 
their trajectories in the phase state are close enough to each other. 

The AI property is understood as a perfect as possible overlap between orbits, histograms and Poincare 
sections (or, ideally, the same - perfect match). 
It can be easy to verify the following property 

Proposition 6 For every N , and Vn, the relation ([s]) can be written in the following form 

N 

Pkirik 

P = • (8) 

k = l 

Moreover, p* is a convex combination of the elements of Vn- 

N N 

The last statement can be easy verified if we denote — m}^/ ^ mj,. Next, because ^ = 1, p* 

k=l k=l 

N 

can be written: p* = ctkPk- 

k=l 

Notation 7 Let denote by A* the synthesized attractor, obtained with the SA implemented by and 
by Ap* the averaged attractor obtained for p = p* . 

Now, the property mentioned at the beginning of this section can be formulated as follows: 
For any N and Ppf, the synthesized attractor A* belongs within the set Aj^. 

The proof is presented in |12l . but the result can be verified by means of computational approach 
too. First, it can be verified computationally that A* and Ap- , with p* given by are AI. Next, using 
Property [6] p* € {pi,pn) , and taking into account the bijection between V and A, we are entitled to 
consider that the same convex structure is preserved from Vn in -An- Then Ap* G An and therefore A* , 
which is AI to Ap-, belongs to {Ai,An) ■ 

Remark 8 i) The time subintervals Ati and the size of the integration step h are parameters which may 
influence the results due to the convergence properties of the considered method for ODEs. Therefore, after 
extensive simulations, we have chosen the best possible values of h so that the best overlap is obtained. 
However, h is not a critical parameter. Therefore, in almost all AS applications, we chose usually values 
for h (in this paper h — 0.005 0.01). 

ii ) To relatively large values for m or N may correspond less or more significant differences between the 
two attractors A* and Ap* ( see 'I I ), but A* remains within of a relatively thin neighborhood of Ap* . 
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The pseudocode of periodic SA, applied on / = [0, Tmax], for chosen A^, Tmax, mi, . . . , mjv, , . . -Pn, is 
presented in Table [l] 

Due to the mentioned convex property, the scheme ([t]) may be applied in any random way |21j . 
Therefore the random SA generates again an attractor A* which, based on the mentioned above convexity 
property, will obviously belong inside the set of considered attractors An endowed with the order of Vn- 
The pseudo-code of one of the possible variants is presented in Table [2] where rand means some random 
generator (in this paper with uniform distribution) of positive integers less than or equal to N. 
m[ count Pi. p* is determined with the following formula 

N 

E Pkm'k 

P* = ^ • (9) 

In this case, in order to obtain a better AI, the integration steps number should be taken as large as 
possible. 

Remark 9 SA cannot be considered as a "true" control algorithm (see Section 3) even it may generate 
any stable trajectory for a considered system, since before the algorithm starts, the system may evolve 
stable and then the algorithm just changes the behavior from a stable attractor to another one. The algo- 
rithm can be use as chaotification algorithm too, but again it should not be considered as a real anticontrol 
algorithm (further informations on chaos control can be obtained e.g. from \2 0j .f23 t/ and for anticontrol of 
discrete and continuous dynamical systems and \25f respectively). The only condition for both control 
and anticontrol is that Vn contain values corresponding to chaotic and stable attractors too. It should be 
notified that the SA can be useful when a desired value for p cannot be set directly. 

To see how the SA must be implement in practice, let us consider the sets Vn and An and suppose that 
certain targeted value of p ^ Vn cannot be accessible, but we want to generate the underlying attractor. 
By using the bifurcation diagram for the considered dynamical system, the only sufRcient condition on p 
is to belong to the real interval {pi,pn) (p cannot be chosen outside this interval because of the mentioned 
convexity property). In order to synthesize the attractor Ap, we must choose rrii so that the desired value 
pis given by Q. This implies to solve ([s]) considered as an equation for fixed Pi, . . . ,Pn, with p* — p, and 
unknowns TOj. With the obtained m.i values, scheme ([t]) can next be applied. The synthesized attractor 
A* will be identical, as shown above, to Ap. Thus, by using the SA, one can "force" the system to evolve 
on the desired trajectory corresponding to p. 

Another practical situation is also possible: Vn and m are not known a priori. Thus, m and the set 
Vni have to be determined so that relation ([S]) be verified with known p. 

In both cases, the solutions are not unique because the elements of Vn belong in a compulsory way to 
one of the infinite number of p-intervals which may compose V . 

For example, let us consider the Lorenz system with the control parameter p, and suppose we want to 
synthesize, with the scheme [toiP2, ^2^1] for fixed h, a stable trajectory corresponding iop= 150 starting 
from Vn — {130, 190}. Then, one of the possible solutions to (|8| is mi = 2 and 7712 = 1 for which ([S]) is 
verified: p* = 150 = (2 * 130 + 1 * 190)/(2 + 1). 
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3 Finding the Rikitake's attractors 

For the Rikitake system modeled by the equations Q, we have 
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Throughout this paper the parameter a is set to the value a = 1. 

In order to apply SA, the bifurcation diagram (Figj2| is a useful tool to study the character of the 
attractors. 

In the same figure, the attractors used in achieving the synthesis are plotted together with the synthe- 
sized attractors. 



3.1 'Transient attractor' 

Numerical simulations of the Rikitake system suggest that the system has attractors which are obviously 
bounded. In other words, the solutions enter a ball around the origin from which they never escape. 
However, this seems not to be true, since if we take as an initial condition (0, 0, k) the exact solution of 
([1]) is 2:i(i) — X2{t) = 0, X3{t) = t + k. But this solution is unstable, the a:3-axis being an invariant 

manifold. There are orbits which escape to, or come from, infinity, instead of going towards the attractor. 
If the flow is on the x^-axis it never escapes; if the flow is not on the x^-axis, then it can never enter. 
Because of the mentioned symmetry of solutions, for points arbitrary close to the x^-axis, the flow takes 
the trajectory back to a bounded attractor (results on boundedness of solutions to third-order nonlinear 
differential equations can be found in |26j). 

Even though for the above mentioned reasons in some papers the x^-axis is regarded as unimportant to 
the dynamics of the system, we found interesting dynamics due to a:3-axis instability. Thus, for p = 90, we 
have found a new and intriguing case which can be considered a kind of "transient attractor" TA (see 
Figj3](a) where the three-dimensional plot is shown and Figj3](d-f) where the phase projections have been 
drawn). TA is actually only a kind of extremely long time transient (its existence being for t < t' with 
t' close to 6 X 10^) before the trajectory reaches the 'real' attractor, a stable limit cycle denoted L (see 
Figjs] (b) where one can sec the limit cycle L whose magnified three-dimensional phase plot is shown in 
Fig. 5 (c)). The dashed lines in Fig[3](f) indicate the well-defined movement sense of TA. The periodic 
characteristics of L can be observed in the three time series, corresponding to xi, xi and X3 depicted in 
Figg (g), (j) and (m) for < e [O, 10^]. The dashed line indicates the moment t = t' , when L is born. 
The details, for < < i' : t e [9.65 x 10^9.75 x 10^] are presented in Figjs] (h), (k) and (n). For 
t > t', the TA transforms into L . Magnified details, D2, for t £ [7 x 10^,7.1 x 10^] are presented in 
Figjs] (i),(l),(o)). Before TA ends in L, it crosses it several times. It can be seen that the size of TA is of 
the order of 10'^ as compared to the small size of the L which is nearly 1000 times smaller. Because a part 
of the trajectory lies on the a:3-axis (or is very close to it (see e.g. Figjs] (d),(e),(f)), the distance being in 
this case not highlighted by the numerical method), we suspect that the TA appears because of a:3-axis 
instability. On this path along the a:3-axis, the speed of the TA is very slow as compared to the loop 
speed (see the vertical peeks in the time series in Figjs] (g),(j),(m) and the details Di in Figjs] (h),(k),(n)). 
Being a stiff system, not all numerical methods have proved to be adequate to obtain the TA. However, 
the many simulations and the TA geometric symmetry, encouraged us to consider TA as being not just 
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some "false" trajectory due to numerical integration, but a real representative component of the insight 
dynamics of the Rikitake system. 

3.2 Attractors synthesis 

By applying the deterministic or random SA, any attractor of Rikitake's system may be synthesized. The 
standard Rungc-Kutta method has been utilized and the most representative cases has been considered. 
Thus, Vn has been chosen so that all kind of behaviors (regular and chaotic) may be considered. 

Using the deterministic scheme [lpi,lp2] with pi — 9.66 and p2 = 12, one obtains the synthesized 
attractor yl* (Figj4] (c)). Figj4] (a) and (b) present Ap^ and Ap^. A* is identical with Ap* for p* given by 
([S]) p* = [pi +P2T/2 = 10.83 (see Figjt] (c) where both A* and Ap. are plotted superimposed and FigjS] 
(d) where histograms of both attractors are presented). It should be noticed that, in this case, SA can be 
viewed as control algorithm (see Remark [oji) since Ap^ and Ap^ are chaotic and A* is a stable limit cycle 
(see FigjI). 

An attractor can be obtained within several variants of ([7|. For example Aio.sa synthesized bellow, can 
be obtained too with the scheme [Ipi, 2p2, Ipa] with = 5, P2 = 7 and p^ = 24.32 (see Fig[2]and Figj5|. 
Here, p* = 10.83 = (pi + 2p2 +P3)/4. A* and Ap* are plotted superimposed in Figj5](d). Superimposed 
histograms (Fig[5](e)) underline the identity. 

With the scheme [Ipi, lp2] with pi — 17 and p2 = 23 the chaotic attractor A* , which is identical to 
Ap, with p* = 20 (Fig(2]and Fig[6] (c)), can be obtained. Superimposed histograms and the Poincare 
section (Figj6](d),(e)) demonstrate the identity. Here SA could be considered as an anticontrol algorithm 
(see Remark]9]i). 

If between pi and P2j corresponding to stable attractors, there are no chaotic windows, anticontrol 
cannot be realized (Remark [9] i) . For example, by choosing pi — 26 and p2 — 28 for which Ap-^^ and Ap^ 
are stable limit cycles (see Figj2]and Figj7](a) where both attractors are plotted together), whatever kind 
of scheme ([t]) is used, only stable limit cycle can be obtained. For example, with the pi.2 chosen bellow, 
and the scheme [lpi,2p2] one obtains the limit cycle depicted in Figj7](b). The identity between A* and 
Ap* with p* = 27 is underlined by histograms and Poincare sections superimposed in FigjTjjd) and (e). 

As stated in Section 2, attractors can be synthesized by using random SA. Thus, by choosing the 
uniform random distribution with algorithm presented in Table |2j for pi — 14.5 and p2 — 20.7, one 
obtains A* identical with Ap* with p* = 17.6 given by ^ (Figj2]and Figjt];c)). 

The results are presented in Table |3] 

4 Conclusions and further directions 

In this paper we verified numerically and computationally that any attractor of the Rikitake system can be 
obtained with the SA, following some deterministic or random rule, the analytical proof being presented 
in [H]. 

In addition to the fact that the SA may force the system to evolve along any attractor, it introduces 
a convex structure inside the attractors set. 

Also, SA may serve as a model to explain the born of specific dynamics in systems encountered in the 
real world or in experiments where this kind of deterministic or random parameter switches may occur. 

In this study we extended the application of the SA from systems with hyperbolic equilibrium to 
systems with non-hyperbolic ones. We have also found an interesting transient attractor, TA, which 
persists for a long period of time. Its presence is supposed to be due to the instability along the x^-axis, 
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and not all numerical methods for ODEs can reveal it. Representing an interesting component of Rikitake's 
dynamics, the TA worth to be studied rigorously, for example from the point of view of shadowing theory. 

The rigorous proof of the existence of the bijection between the set of parameter values and the set of 
attractors and the analytical proof of the application of the SA to other class of systems represent tasks 
for future studies. 
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Table 1: Pseudo-code of the PS algorithm 



repeat 

for i = 1 to mi do 



end 



integrate ([2j) for p = 

t = t + h 



for i = 1 to m^v do 

integrate ^ for p = 

t = t + h 

end 

until t > Tmax 



Table 2: Pseudo-code of the random PS algorithm 

repeat 

label=rand(iV) 
if label = 1 then 

integrate ^ with p = p^ 

nil — ™i + 1 
if label — 2 then 

integrate ^ with p = p^ 

m'n — m'o + 1 



if label ~ N then 

integrate ^ with p = pn 
m'jy = + 1 
t = t + h 
until t > Tmax 
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Table 3: The results of PS algoritm applied to the Rikitake system ([T]) 



Scheme 


PI 


P2 


P3 


P* 


Remarks 


[lpi,lP2] 


9.66 


12 




10.83 


Apj , chaotic, A* stable hmit cycle (Fig 


4 




[1P1,2P2,,1P3] 


5 


7 


24.32 


10.83 


Apj , Ap2 chaotic, Ap^ stable limit cycle, A* stable limit 


c 


ycle (Fig 5 I 


[1P1,1P2] 


17 


23 




20 


Apj^ , Ap2 stable limit cycles. A* chaotic (Fig 


6 


1 


[1P1,1P2] 


26 


28 




27 


Apj , Ap2 stable limit cycles, A* stable limit cycle (Fig 


7 


) 


random SA 


14.5 


20.7 




17.6 


Apj , Ap2 chaotic attractors, A* stable limit cycle (Fig 


8 





12 




13 




Aal A» 




*M! MM !tli!M!ni !!! ! i M ! 1 ! ! ! ! Ml | ■ ! I ■ ! ! 




(0 



Figure 2: Bifurcation diagram for the Rikitake system, (a) Bifurcation diagram for the component 
Xi where the positions of the attractors utihzed and the synthesized attractors are indicated; (b), (c) 
Bifurcation diagram for components X2 and X3. 
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Figure 3: Transient attractor (TA) and stable limit cycle L. (a) Phase plot; (b) detail; (c) Limit cycle L 
after neglected TA; (d), (e), (f) projections of TA and L; (g) Time series for xi component; (h), (i) the 
details Di and D2 indicated in (g); (j) Time series for component X2; (k),(l) details Di and D2 indicated 
in (j); (m) Time series for component x^; (n),(o) Details Di and D2 indicated in (m). 
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Figure 4: Synthesized attractor A* and Ap* obtained with scheme [Ipi, IP2] with pi = 9.66 and p2 = 
12; p* = 10.83. (a), (b) Phase plots of Ap^ and Ap^- (c) A* and Ap. superimposed; (d) Histograms 
superimposed. 
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Figure 5: Synthesized attractor A* and Ap. obtained with scheme [Ipi, 2p2, Ips] with pi = 5, p2 = 
7, P3 = 24.32; p* = 10.83. (a)-(c) Phase plots of Ap-^ 3 (d) A* and Ap. superimposed; (e) Histograms 
superimposed. 
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Figure 6: Synthesized attractor A* and Ap* obtained with scheme [lpi,lp2] with pi — 17, p2 — 23; 
p* — 20. (a), (b) Phase plots oiAp^ and Ap^; (c) A* and Ap. superimposed; (d) Histograms superimposed; 
(e) Poincare sections superimposed. 
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x1 range 

(c) (a) 

Figure 7: Synthesized attractor ^* and Ap* obtained with scheme [lpi,lp2] with pi — 26, p2 — 28; 
p* = 27. (a) Phase plots of Ap-^ and Ap^ ; (b) A* and Ap* superimposed; (c) Histograms superimposed; 
(d) Poincare sections superimposed. 
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Figure 8: Synthesized attractor A* and Ap* obtained with random scheme with pi = 14.5, p2 = 20.7; 
p* = 17.6. (a), (b) Phase plots of Ap^ and Ap^; (c) A* and Ap-* superimposed; (d) Poincare sections 
superimposed; (e) Histograms superimposed. 
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